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Abstract 

We investigate the thermodynamic Bethe ansatz (TBA) equations for a 

system of particles which dynamically interacts via the scattering matrix 

of affine Toda field theory and whose statistical interaction is of a general 

i) , Haldane type. Up to the first leading order, we provide general approximated 



"^ ' analytical expressions for the solutions of these equations from which we 

derive general formulae for the ultraviolet scaling functions for theories in 

^ . which the underlying Lie algebra is simply laced. For several explicit models 

H I we compare the quality of the approximated analytical solutions against the 

numerical solutions. We address the question of existence and uniqueness 
of the solutions of the TBA-equations, derive precise error estimates and 
determine the rate of convergence for the applied numerical procedure. A 
general expression for the Fourier transformed kernels of the TBA-equations 
allows to derive the related Y-systems and a reformulation of the equations 
into a universal form. 
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1 Introduction 

More than twenty years ago Yang and Yang [||] introduced in a series of seminal 
papers a technique which allows to compute thermodynamic quantities for a system 
of bosons interacting dynamically via factorizable scattering. This method was 
generalized less than ten years ago by ALB. Zamolodchikov Q] to a system of 
particles which interact dynamically in a relativistic manner through a scattering 
matrix which belongs to an integrable quantum field theory. This latter approach, 
usually referred to as thermodynamic Bethe ansatz (TBA), has triggered numerous 
further investigations [3-13]. The reason for this activity is twofold, on one hand 
the TBA serves as an interface between conformally invariant theories and massive 
integrable deformations of them and on the other hand it serves as a complementary 
approach to other methods. 

The TBA-approach allows to extract various types of informations from a mas- 
sive integrable quantum field theory once its scattering matrix is known. Most 
easily one obtains the central charge c of the Virasoro algebra of the underlying 
ultraviolet conformal field theory, the conformal dimension A and the factor of 
proportionality of the perturbing operator, the vacuum expectation values of the 
energy-momentum tensor {T^\ and other interesting quantities. Thus, the TBA 
provides a test laboratory in which certain conjectured scattering matrices may be 
probed for consistency. 

In addition, the TBA is useful since it provides quantities which may be em- 
ployed in other contexts, like the computation of correlation functions. For instance 
the constant of proportionality, the dimension of the perturbing field and (T'^) may 
be used in a perturbative approach around the operator product expansion of a two 



point function within a conformal field theory ||I^. (T'^) may also be used as an 
initial value for the recursive system between different n-particle form factors |]n . 

In order to obtain data beyond c one has to investigate the behaviour of the 
scaling function c{r) (r is the inverse temperature times a mass scale), which may 
be viewed as the deformed central charge of the Virasoro algebra. It will be the 
central aim of this manuscript to determine this function. Certain statements in 
this context can be made in complete generality, but of course ultimately one has 
to specify some concrete models in order to be more explicit. An attractive choice 
for these models are afiine Toda field theories [jT5[, since they cover a huge class 
of theories due to their Lie algebraic formulation and permit therefore to extract 
many universal features. Ultimately they also allow a comparison against standard 
perturbation theory in the spirit of [^ . 

The first computation of such a scaling function was carried out in for the 
scaling Lee- Yang model and scaling 3-state Potts model (minimal A2 -affine Toda 
field theory). Meanwhile there exist several computations of this kind f^, H, 0, p!T|| . 



The general behaviour one finds is 



c{r) = Ceff + /(r) + Y: c„r2"(^-^) 



where c^s = c — 24h' is the effective central charge, with c being the usual conformal 
anomaly, i.e. the central charge of the Virasoro algebra, and h' the lowest conformal 
dimension of the underlying conformal field theory. The c„ are some constants 
and typically /(r) ~ r^ or /(r) ~ r^ln(r). However, in |jl2|, |l^ a quite different 



behaviour was observed, namely that /(r) ~ (const — ln(r))^^, which was attributed 
therein to zero mode fluctuations. It is this kind of behaviour which we find for all 
affine Toda field theories (see equn. (^)). We like to stress that the constant const 
should be different from zero. 

In our investigation we slightly modify and generalize an approach which was 



introduced originally in [12|. We also include in our analysis the possibility that 



the statistical interaction is of general Haldane type |0] and investigate the TBA- 



equations adapted to this situation [|1^]. We investigate in detail the expressions 
for the TBA-kernels and their Fourier transformed versions, from which we obtain 
a universal form for the relevant TBA-equations. We derive the related Y-systems. 
Furthermore, we address the question of the existence and uniqueness of the solu- 
tions of the TEA equations, derive error estimates and the rate of convergence for 
the numerical procedure applied. Our analytical considerations culminate with the 
derivation of a general expression for the scaling function valid for all affine Toda 



field theories in which the underlying Lie algebra is simply laced (|83|) and which 
depends in a universal manner only on the rank of the algebra, its Coxeter number 
and the effective coupling constant. We compare the general formulae against the 
numerical solutions for several explicit models, i.e. for the 74^,742,^42,^43 and 
{G2 , -D4 )-affine Toda field theories. 

Our manuscript is organized as follows: In section 2 we recall the TBA-equation 
for general Haldane statistics and explain how it can be used to compute explicitly 
expressions for scaling functions. In section 3 we derive general expression for the 
approximated analytical expression for the solutions of the TBA-equations and for 
the scaling functions c{r). In section 4 we concretely test our statements for a 
system involving the scattering matrix of affine Toda field theories. We compare 
the general approximated analytical expression with the numerical solution for some 
explicit models. In section 5 we address the question of the existence and uniqueness 
of the solutions of the TBA equations and derive error estimates and the rate of 
convergence. We state our conclusions in section 6. 



2 The TBA-equations 

We consider a multi-particle system involving / different types of particles and 
assume that the two-particle scattering matrix Sij{6) (as a function of the rapidity 
6) together with the corresponding mass spectrum {rrii is the mass of particle type 
i) have been determined. The / coupled TBA-equations which describe such a 
system, based on the assumption that the dynamical interaction is characterized 
by the scattering matrices Sij{d) and whose statistical interaction is governed by 
general Haldane statistics [^ (^jjQ, read |Tr 



cosh^ + In (l - e-^"(^'")) = J2 ("^u * Lj) {6, r) . (2) 



The scaling parameter r is the inverse temperature times a mass scale and we denote 
as usual the convolution of two functions / and g hj {f * g){0) := l/(27r) / dO' f{6 — 
6')g{6') and the kernel by 

^M := -z-^ In S,,{9) - 2ng,,6 (9) = ^,^{9) - 2ng,,6 (9) . (3) 

We recall the important fact that the functions Li{9,r) are related to the ratio of 
densities p^of particles inside the system over the densities of p\ of available states 
as Li{9,r) = ln(l + pi./ pi). This implies that Li{9,r) > 0, a property we will 
frequently appeal to. In addition we will make use of the fact that Li{9,r) is an 
even function in 9. Once the system of coupled non-linear integral equations (0) is 
solved (usually this may only be achieved numerically) for the / unknown functions 
Li{9,r), one is in a position to determine the scaling function 

or . r 

cir) = — Y,m J d9L,{9,r) cosh9 , (4) 

i=l 

which may be viewed as the off-critical central charge of the Virasoro algebra. We 
shall now adopt the method of [|12| and instead of regarding the TBA-equation as 
an integral equation, we transform it into an infinite order differential equation. 
This is possible with the sole assumption that the Fourier transform of the function 
^ij{9) can be expanded as a power series 

<l>,j{k) = d9<^ij{9)e'''^ = 27c Y.{-i)''r]\fk'' . (5) 

jL "=0 



*The bosonic and fermionic statistics correspond to gij — and gtj — Sij, respectively. 



Then it is a simple consequence of the convolution theorem^ that the integral equa- 
tion (^ may also be written as an infinite order differential equation 

I oo 

rm, cosh + In (l - g-^'^^-'-)) = ^ E ^1^4"^' ^) • (6) 

j=l n=0 

We introduced here the abbreviation Lf" \0^r) = {d/d9)'^Li(9,r). The two alter- 
native formulations of the TBA-equations serve different purposes. Its variant in 
form of an integral equation (|]) is most convenient for numerical studies, whereas 
the differential equations (^ turn out to be useful for analytical considerations. 

From a numerical point of view it is advantageous to make a change of variables 
and instead of the functions Li{9, r) use the so-called pseudo-energies ej(6', r) defined 
via the equations 

e.(^,r) :=-E^.,L,(0,r)-ln(l-e-^»(^''^)) . (7) 



For fermionic and bosonic type of statistics the formulation of the TBA-equations 
in terms of pseudo-energies is most common. However, since (0) may obviously not 
be solved in general, it is more convenient to keep the TBA-equations in terms of 
the Li[6,r) for general Haldane type of statistics. Alternatively we may relate all 
types of statistics to the fermionic one in the following way. Considering a system 
which interacts statistically via gij, we can define the quantity g'-- = 6ij — Qij and 

parameterize all L— functions by Li(9,r) =: Ci{9,r) = Infl + e"^**-^''"'' j. Then the 
TBA-equations (0) may be rewritten as 

6,(9, r) = rm, cosh 9 - ^ [((y.,, + 27ig:^S) * £,) {9, r)] . (8) 

i=i 

It should be kept in mind that 6i{9,r) {^ ei{9,r)\) is now a formal parameter 
and, except in the fermionic case, it is not related to the ratio of particle densities 
pl/Ph in the characteristic way as a distribution function associated to the relevant 
statistics, i.e. ei{9,r) obtained as a solution of the equations (0). 

3 The ultraviolet Limit 

We shall now analytically investigate the behaviour of the scaling function in the 
ultraviolet limit, i.e. r is going to zero. For this purpose we generalize the procedure 



of Zamolodchikov |T2|, which leads to approximated analytical expressions. We 



attempt to keep the discussion model independent and free of a particular choice 



of the statistical interaction as long as possible. Introducing the quantity Li{0) : = 
Li{6 — ln(r/2), r) and performing the shift 9^9 — ln(r/2), the TBA-equations (P) 
acquire the form 

I oo 

m,e' + In (l - e-^'(^)) = ^ E vfLf\9) . (9) 

j=l n=0 

Here we have neglected the terms proportional to e^''^^'"/^)"^, under the assumption 
that 21n(r/2) -C 9. Obviously the r-dependence has vanished, such that the Li{9) 
are r-independent. The equation for the scaling function @ becomes, under similar 
manipulations and the neglect of similar terms as in the derivation of (j^) , 

I 



6 

) 

" 4 = 1 



cW = ^E^^ / d9U[9)e' . (10) 



ln(r/2) 



In analogy to the procedure of |T^, we consider now the so-called "truncated 
scaling function" 

c{r,r') = —J2m^ld9Li{9)e' , (11) 

which obviously coincides with c{r) for r' = ln(r/2). At this point we make several 
assumptions: 

i) The functions Li{9) obey the equations (P). 

ii) In the power series expansion (^) the coefficients are symmetric in the particle 
type indices, i.e. rjlj = rij^ . 

iii) All odd coefficients vanish in (|^), i.e. rjlj"' = 0. 

iv) The asymptotic behaviour of the function Li{9) and its derivatives read 

lim Lf V) = forn>l,l<i</ , (12) 

lim e'^Li{9) = for 1 < i < / . (13) 

Assumption ii) is guaranteed when the two-particle scattering matrix is parity 
invariant and the statistical interaction is symmetric in particle type, Qij = Qji. 
The requirement iii) puts of course constraints on the scattering matrices for given 
statistics, or vice versa. It will turn out in the next section, that for fermionic 
statistics it is satisfied by all scattering matrices of interest to us. The asymptotic 
behaviour iv) will be verified in retrospective, that is all known numerical solutions 
exhibit this kind of asymptotic behaviour. 



Under these assumptions the truncated scahng function may also be written as 

q / / CO 271—1 \ 

^ i,j=l \n=l k=l I 

-^ E (^(1 - e-^^^^'^) + ^ ln(l - e-^^('^')) + m^e^'uA . (14) 



Here L{x) = J2'^=i x"'/n'^ + l/2 ln(a;) ln(l — x) denotes Rogers dilogarithm (e.g. [pJ| ). 
The equahty ( [T^ ) is most easily derived when considering the differential equation 
which results from (pA]) 

^l^.-l,e"±,nMr') • (15) 

2 = 1 

One may now verify by direct substitution^ that (|15D is solved by (|l^) when i)-iii) 



hold. The constant of integration is fixed by the property c{r, oo) = 0, such that 
(|I^) is exact if iv) holds. 

3.1 The extreme Limit 

One of the best known outcomes of the TBA-analysis is the fact that in the extreme 
ultraviolet limit the scaling function becomes the effective central charge of the 
underlying ultraviolet conformal field theory, i.e. lim^^o c{f^) = Cea- In order to 
carry out this limit, we note that the assumptions (|T^ and (|T^) also imply 

lim Lf ^(^, r) = forn > 1 , 1< i < / , (16) 

such that equations (|^) become a set of coupled non-hnear equations for the / 
constants Lj(0, 0) 

ln(l-e-^'(°'°))=Ei°^^.(0,0) . (17) 

Due to the fact that Li{9,r) is positive, we deduce that these equations admit 
physical solutions when 7]}^ < 0. For a given S-matrix this condition will restrict 



possible statistical interactions. With the help of the solutions of (|T^ we recover 
from ([I^) the well known formula for the extreme ultraviolet limit 

limc(r)=Ceff=-^EL(l-e-^»(°'°)) . (18) 

2=1 

These equations have been analyzed extensively in the literature for various models 
0, 1^, ^, ll^] and as we demonstrated, also hold for general (Haldane) type of 



statistics. 



■I-The identity / dt ln(l — e *) = i(l — e ^) + x/2 ln(l — e ^) is useful in this context. 





3.2 Next leading Order 

In order to keep the next leading order in the ultraviolet limit we proceed as follows: 
Under the assumption that the symmetry property of Li{6) still holds at this point of 
the derivation, we may neglect in (|^) also the term e^. Assuming that Li{6) is large 
and ([T2I) holds, the TBA-equation in the ultraviolet limit may be approximated by 

E(r/lf^fW+r/l?4W) + e-^'(^) = . (19) 

Unfortunately this equation may not be solved analytically in its full generality. 
However, choosing now the statistics in such a way that the rjl = Z]j=i Vij = we 
may solve (|19|). In most cases we shall be considering in the following, this condition 
implies that we have fermionic type of statistics. For this case the solution of ([19|) 
reads 

with Oj, /5j, Oj, /3j being the constants of integration and rji = J2j=i Vij ■ Note that 

(2) 
the assumption ^j 7^ is not always guaranteed below. We will discard the second 

term in the following w.l.g. Under the same assumptions the expression for the 



truncated scaling function (]1^) becomes 

3 ' 



c(r,r') = / + 4E(^r(4^'HO)'-2e-^'(^')) . (21) 

Substitution of the solution (|^) into (^l]) yields 

e(r,r')=/-^Er^fV . (22) 



(2)^,2 
1=1 



Notice that this expression for the truncated effective central charge is independent 
of the constants f3^ and r'. We use the latter property to argue that in fact the r.h.s. 
of (p2D corresponds to the scaling function c{r). Invoking now also the property 
Li{6) = Lj(21n(r/2) — 6) we obtain the additional relations 

nvr 
"' = m - ln(./2)) • <'"' 

with n being an odd integer, which we choose to be one. At the moment we do not 
have any further argument at hand in order to fix the remaining constant. Therefore 
we have 

He,r) ^ ..(^^) (24) 



Lf\e,r) = -2aitan(M) (25) 

L^^\e,r) = -2a^J cos^ (aiO) (26) 

Lf{9,r) = -Aaftan{a,9)/cos^{aie) (27) 

Li-\9,r) ~ < . (28) 

Using the fact that a^ tends to zero for small r, the equations (|25|)-(|28|) demonstrate 
the consistency with the assumption that the derivatives of Li{6, r) with respect to 



6 are negligible, i.e. equation (p!6|). Closer inspection shows that for given r the 
series build from the Lf {O^r) starts to diverge at a certain value of n. Since (|2^) 



is not exact this does not pose any problem, but one should be aware of it. The 
scaling function becomes in this approximation 

~{ Wi - ln(r/2))2 

As already pointed out there does not seem to be any argument in this approach 
which allows to fix the constant (3^. However, we will present below a natural guess 
and also resort to numerical data to fix it. 

In order to perform more concrete computations one has to specify a particular 
model at this point. 



4 Affine Toda Field Theory 



Affine Toda field theories [T^ constitute a huge, important and well studied class 
of relativistically invariant integrable models in 1+1 dimensions. The exact two- 
particle scattering matrices of all affine Toda field theories with real coupling con- 
stant were constructed on the base of the bootstrap principle [^ and its generalized 



version |2^. Various cases have been checked perturbatively. The theories exhibit 
an entirely different behaviour depending on whether they are self-dual or not. Here 
duality has a double meaning, on one hand it refers to the invariance of the algebra 
under the interchange of roots and co- roots (i.e. A'-^\ A2J, D^\ Eq',E-}',Eq') and 
on the other hand it refers to the strong- weak duality in the coupling constant. 

4.1 The S-matrix 

4.1.1 Simply laced Lie algebras 

We recall now some well known facts about the two particle scattering matrix and 
present some new features which will be important for our analysis below. For 
theories related to simply laced Lie algebras the two particle S-matrices [^, ^, ^ 



involving the r different types of particles may be furnished into the universal form 

S^A0) = X{\2q-^^-^^\ , l<z,j<r = rankg. (30) 

g=l I ^ ie 

The building blocks {x}q may be expressed in terms of hyperbolic functions, albeit 
it will be most convenient for our purposes to use their integral representation 

^ ^ 



with 



/.,B(t) = 8 sinh ('^) sinh (^ (2 - B)\ sinh {^j^{h-x)] . (32) 



We adopt here the notations of [^ and denote the Coxeter number by h, funda- 
mental weights by Aj, the colour values related to the bicolouration of the Dynkin 
diagram by c{i) and a simple root a^ times c{i) by 7^. The Coxeter element is chosen 
to be 0" = o"-cr+, where the elements a±. = HjeAi (^i are introduced, with ctj being 
a Weyl reflections and A± the set of simple roots with colour values c{i) = ±1, 
respectively. The effective coupling constant B{P) = (/?^/27r)/(l + /3^/47r) depends 
monotonically on the coupling constant (3 and takes values between and 2. 



Many of the properties of ( pOD are very well documented in the literature [0, 
|22| , |2^ and we will therefore only concentrate on those which are relevant for our 
investigations. One of the remarkable features is the strong-weak duality in the 
coupling constant /3, which is a particular example for one of the dualities which are 
currently ubiquitous in string theoretical investigations. Whilst in the latter context 
duality serves as a powerful principle, we will frequently employ it as a simple 
consistency check on various expressions which arise during our computations. The 
strong-weak duality manifests itself by the fact that each individual building block, 
and consequently the whole S-matrix, possesses the symmetry B ^ 2 — B. A 
further check, which we wish to pursue from time to time for consistency reasons, 
is taking the coupling constant to zero, i.e. -B(O) =0 or equivalently B{oo) = 2, 
such that we obtain a free theory with Sij{6) = 1. Also we want to compare with 
some known results for the so-called minimal part| of the scattering matrix (pO]), 
which is formally obtained by taking the limit B —>■ zocpj. 

It will be important below to recall that the block {l}g only occurs (with power 
one) in the scattering matrix between two particles of the same type, i.e. in Sii{d), 



§ This part satisfies by itself the bootstrap equations and the additional factor is of a CDD- 
natiirc, meaning that is does not produce any relevant poles inside the physical sheet. 

^ For the blocks {x}g in form of hyperbolic functions this limit is trivial. For the integral 
representation (JS^ ) it requires a bit more effort, but is easily verified with the help of the Riemann- 
Lebesgue theorem (If 17(3:;) e ii(—cxD,cx3) then lim /_ g{x)e~'^*^dx — 0). 



which imphes for the S-matrices (pO| ) Sij{0) = (— 1)"^'^. A further relation which 
exhibits the occurrence of particular blocks is 

S., (0) = {2}i^ n 2^? - ^ , ^^^ , c{i} + c{j) . (33) 

q=2 { '^ ) 

Here lij is the incidence matrix of the Lie algebraic Dynkin diagram, i.e. twice the 
unit matrix minus the Cartan matrix K. Since relation (|33| ) will be important for 
our analysis and does not seem to appear in the literature, we will briefly prove it 
here. We make use of the interrelation between the incidence matrix and particular 
elements of the Weyl group and the action of these elements on simple roots El, EH] 



Y, hj>^j = {(T- + CT+)Ai and ae^a^ = -a^^^^'^-'^'^'^^^aj . (34) 

i=i 

Taking the inner product of the first equation in (l3^) with aj and the subsequent 
application of the second relation on the r.h.s. yields 

lij = Aj ■ cT-c{j)<yj + Aj ■ (Tc{j)aj = -Aj ■ aj - Aj ■ a'^'-^^'aj . (35) 

Therefore, by the orthogonality of fundamental weights and simple roots and the 
symmetry in i and j, it follows directly that 

A,-a% = T/,,, c{{}^c{j) . (36) 

Noting further that {x}g = {2h — x}g , we may extract the block {2}^ with power 
lij from the product in ( ppj ) and obtain (|5B|). 

Remarkably one may also carry out the product in ( pOD and obtain a closed 
expression for 5* in the form 



(oo , 



/(2i)),B(t)(2cosh^-/y -2%jsin(^^)| . (37) 



The Lie algebraic quantities involved in this identity, i.e. / and /i, are more easily 
accessible than the orbits of the Coxeter elements, however, this is at the cost that 
the singularity structure is less transparent. Formula ( P7| ) coincides with equation 
(5.2) in p5[ (up to a factor (—1)"^*^), once the general expression in there is reduced 



to the simply laced case. 

In the course of our argumentation we will also employ the identity 



s,, ( ^ - y ) s.,, [e + '^^ = \{ Su (ef^ MO, (38) 



which was derived first in 0. 

10 



4.1.2 Non-simply laced Lie algebras 

Once one turns to the consideration of afiine Toda field theories related to non- 
simply laced Lie algebras many of the universal features which one observes for 
simply laced Lie algebras cease to be valid. For instance for the simply laced Lie 
algebras it is known that the singularities of (|30|) inside the physical sheet do not 
depend on the coupling constant and remarkably all masses renormalise with an 
overall factor. These latter behaviours dramatically change once the related Lie 
algebras are taken to be non-simply laced (except Ag^ which is also self-dual). 
In fact in these cases one can not associate anymore a unique Lie algebra to the 
quantum field theory, but a dual pair related to each other by the interchange of 
roots and co-roots. General formulae similar to (pOD, which are valid for all non- 
simply Lie laced algebras have not been constructed yet. However, one can still 



construct S-matrices 26, 27 



S:,{e) = m^Me (39) 



x,y 



out of some universal building blocks which generalize (|31|) 



^ ^ 

In distinction to the simply laced case, the variable x may now become non-integer, 
y is 0, 1/2, 1 or -1/4 and 

g,,y,B{t) = 8sinh (^(1 - 2y)) sinh (^ (2 - 5)) sinh (|: (H - x)) . (41) 

The Coxeter number h has been replaced by a "fioating Coxeter number " H, in the 
sense that it equals the Coxeter number of one of the two dual algebras in the weak 
limit and the other in the strong limit. The singularities of (|39|) now depend on 
the coupling constant which leads to a modification of the bootstrap as explained 
in [Q. The shift of the masses resulting from the renormalization procedure can 
not be compensated anymore by an overall factor. 

In contrast to the simply laced case, the control over the set in which x and 
y take their values is slightly different than in (^). Whereas in (^) the powers 
of the blocks may be computed directly from the Lie algebraic quantities, in the 
non-simply laced case they are obtained solving at first some recursive equations 
and thereafter extracting the powers from a generating function |2^. There exists 
however a generalization of (|37|) which also includes the non-simply laced case. In 
this paper we will not treat the completely generic case and shall be content with 
treatment of one particular case in order to exhibit the difference towards the simply 
laced case, such that (|39| ) and (pO[) are sufficient for our purposes. 

11 



4.2 The TBA-kernels 

Given the S-matrices we are now in a position to compute the relevant quantities for 
our analysis, i.e. the kernels appearing in the TBA-equations (|^) and the Fourier 
coefficients in equation (|^). 

4.2.1 Simply laced Lie algebras 

Taking the logarithmic derivative of the phase of the scattering matrix (pOD the 
TBA-kernels (without the statistics factor) are easily computed to 

1 ^ 

^iji^) = -9 E(^^ ■ ^''7i)^2g-(c«+c0))/2(^) , (42) 

^ q=l 

where 

u^M = -'^^''i'^}o = lJ^tUB{t)cos(^A, (43) 

^ ^ 

= 7.-i(^)+7.+i(^)-7.+B-i(^)-7.-B+iW- (44) 

We introduced the function 73,(6') = sin(7rx//i)/(cos(7r2;//i) — cosh^) and set 7o(^) = 
0. Depending on the context either the variant ( ^31) or (^) turn out to be more 
convenient. The Fourier transform is most easily evaluated when we use the integral 
representation (^31) for each block 



= lx-ii.k) + l^+Ak) - 7^.+B_i(fc) - 7x-B+i(^), (46) 

where ^^i^) = 27rsinh[(x//i — l)7rA;]/sinh(7rA;). Therefore we have 

^^Ak) = -lj:iK ■ a''+W^-)-^)/27,)c:'2,+(cO-)-eW)/2-i(^) • (47) 

This means, that just like the S-matrix itself, we may express the TBA-kernels 
and their Fourier transformed versions in terms of general blocks. The difference 
is that now instead of the powers of the blocks their pre-factors are determined 
by orbits of the Coxeter element. Notice also that the property {x}g = {2h + x}g 
guaranteed that in ( pO]) the products Ylg=i and YlgZo ^^^ equivalent. Since now 
7a; (^) 7^ l2h+xik), wc have to take the sum over the appropriate range. This is the 
reason why q is shifted in (^^. Remarkably for the Fourier transformed kernels it 



is possible to evaluate the whole sum over q in (|47|) in an indirect way. We obtain 
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vr \ ^ / ^ , TT 



the universal form 

ipij{k) = 27r^('/-2cosh^A;j f / - 2cosh ^A;(l - 5) j (48) 

, f(B-2)nk\ , fBnk\ ( , vrfc \"^ . . x 

= Stt sinh ^ — ^ — sinh — — 2 cosh -^ - ^ - '^T^^ij ■ (49) 

A similar matrix identity has turned out to be extremely useful in the inves- 
tigation H] of the TBA- kernel for the minimal part of the scattering matrix (pOl). 



Equation (^Sj) coincides with an identity quoted in |13[ once the roaming param- 
eter therein is chosen in such the way that the resonance scattering matrices take 
on the form of (|30|) . Taking the inverse Fourier transform of (|49| ) we obtain after 



an intgeration with respect to 9 the integral representation (pTf ) for the scattering 
matrix. 

Since the identity (^) is by no means obvious, in particular with regard to 
(^), we will now provide a rigourous proof of it. Essentially we have to use the 
Fourier transformed version of the relation ( |55| ) for this purpose. Shifting the Fourier 
integral into the complex plane yields 

\im <f deipijie)e"'' = ipijik)-p JdeifiJe±t^)e"'^'^'''/''^ (so) 






= 27T6ije^^^''/'' - Tilije^^^'^ . (51) 

Here the contours C^ are depicted in figure 0. V denotes the Cauchy principal value. 
The poles inside the contours are collected by considering (^) and (PI). Obviously 
coupling constant dependent poles may only result from UiiO) and (^2/1-1 (^), whilst 
poles directly on the contours at ^ = iivr/Zi may only originate from uj2{0) and 
uj2h-2{Q)- From the statements made in subsection 4.1.1., it follows that the former 
blocks may only occur when i = j and from equation (|33|) we infer that the pre-factor 
of Ci;2(^) is lij. The relevant residues are computed easily and by noting further that 
the singularities directly on the contour count half, we have established (|5T|) . 

Acting now with —i times the logarithmic derivative on the identity (PHf), mul- 
tiplying with exp{ik9) and integrating thereafter with respect to 9 we obtain 

V 

V jde (y,,^. {9 + 17, /h) + ^,^. {9 - iT^/h)) e^'' = Y: hm,{k) • (52) 

On the other hand the l.h.s. of ( |52D may be computed alternatively from the right 
hand sides of (|50|) and fPD, such that we obtain 



y^ IiiiPij{k) = 2ip^j{k) cosh —k + 27clij — ArcSij cosh 7-^(1 ~ B) , (53) 

1=1 "- "- 
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Jm{9) 



Tl/h 



ttB 
h 






ttB 
h 



^^FJh 



C: 



Re{9) 



Figure 1: The contours C^ in tiie complex 6'-plane. The bullets • belong to poles resulting 
from UJ2{9) and the open circles o to poles of uJi{0). 



and therefore (0). 

When we take the limit B — > zoo, the coupling constant dependent poles move 
outside the contours C^, such that the term involving B will be absent in (^) and 
we recover the result of [Q. The weak limit B ^ turns equ. (53) into a trivial 
identity by noting that with (^) lim^^o^ijl^) = 27r(5jj. 

As the last quantity which will be important for our considerations we have to 
compute the power series expansion of (f>^Ak). For this purpose we make use of the 
identity 

sinh(a;t) ^ 2^"+^ ^ (l + x\,^ 



sinht 



which follows directly from the generating function for the Bernoulli polynomials 
Bn{x) of degree n (e.g. [^). We may then expand each building block of the 
Fourier transformed TBA-kernel (E5|) as 



uM = 2vr^/ii2")fc2" 



n=0 

oo 



2vrE(ei + 



{2n) (2n) (2n) 



,(2n) 



'^x-i ^x+l-B ^x-l+B I ^ 



n=0 



2n 



(55) 
(56) 



where the coefficients are I'f'^'' = 2(27r)^"/(2n + l)\B2n+i{x/{2h)). Therefore we 
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obtain for the Fourier coefficients in equation (|^) 

/ 1 \n+l h 

JO) _ /).._„.. and 77^^") - ^~^ VfA ,^9+K?Vi)/2.- ^,,(2") /r7^ 

Vjj — Oij y«j ana r^^^- — ^^i^a^ ct 7jJA''2g+{c{i)-c(i))/2-i l-^'J 

^ 9=1 

for n = 1,2,... The second order coefficient is particularly important with regard to 
our approximated analytical solution presented in section 3.2 and we will therefore 
analyze it a bit further. From ( p6D we obtain fi^^^ = n'^B^B — 2){h — x)/h?, such 
that 



h 
(2) 

9=1 



.. - ^ E(A. ■ a'^^^<^^-'^"iyB{B -2){h-2q- {c{j) - c(.))/2 + 1) . (58) 



We may evaluate the sum over q by noting further that 

^^ = r E q^'-'^'^'^-'^^H and fl ^^ " ^'^j = • (59) 

q=l q=l 

The ffist identity follows by inverting the relation 7^ = (1 — a~^)a^^~^^'^^^^'^\i [^ 
and the second by computing the geometric series. Therefore ( |58D becomes 

(2) _ fB{2-Bl _ n^B{2-B) 

^ij ~ /^2 ^» ^J - p ^ij ■ ^^^> 

The latter formula follows immediately by recalling that Aj = I]j=i -^rj^^j- It has 
the virtue that it is directly applicable since it involves only quantities which may 
be effortlessly extracted from Lie algebraic tables. At last we may compute the 
quantity ri\ = Yl^j^iVij which occurs in our approximated analytical approach 
(pop. Obviously we obtain from (pO| ) 

where p = J2l=iK = 1/2Z!q>oQ^ is the Weyl vector. The inner product of Aj ■ 
p may be related to a universal quantity, namely the index of the fundamental 
representation Aj 

dimAj , , - , , s 

^A, = ^^(2p + A.)-A., (62) 

such that we finally obtain 

(2) 
Hence we are also able to evaluate rji in a universal manner. We may proceed 

further and also compute 

^ (2) ^ Tr^B{2-B) 2 ^ Tr^B{2 - B)Q^dimg ^ n'B{2 - B)r{h+ 1) 

^''^ h^ ^ 24/i2 I2h ■ ^ ' 
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Here we used the Freudenthal-deVries strange formula p^ = Q^dim.g/24: (see e.g. 
29[| ) , the fact that for simply laced algebras the eigenvalue of the quadratic Casimir 



operator is Q^ = 2h and that we have dimg= r{h + 1) for the dimension of the Lie 
algebra g. 

4.2.2 Non-simply laced Lie algebras 

For the non-simply laced cases we may proceed similarly, just now we are lacking the 
universality of the previous subsection. Using the generalized blocks (|40|) instead 
of (|3TD we derive 

^ ^ 

= 7x^yB-l{0) + ^x+yB+1 (^) " lx+yB+B-1 (^) " lx-yB-B+1 (^) • (66) 

The Fourier transformed of ujx^y{9) reads 

~ ,,. -^ 9x,y,B{Trk) 

smh(7rA;) 

oo 

= 27rE/^S^fc'" (68) 

n=0 

oo 

— ^^ Z^ \'^x+l+yB ^ '^x-l-yB '^x+l-B-yB '^x-l+B+yBJ '*' ' K^-') 



n=0 



In particular 



, (2) _ Bn'jB - 2)iH - x)il - 2y) 



is important for our purposes. 

4.3 The TBA-equations 

Having computed the TBA-kernels, it appears at first sight that, apart from a 
simple substitution, there is not much more to be said about the form of the TBA- 
equations. However, Zamolodchikov observed in |^ that once the Fourier trans- 
formed TBA-kernels admit a certain representation (in our situation this is (|48|)), 
the TBA-equations may be cast into a very universal form. 

We will exploit now the universal features derived in the preceding section. 
In order to keep the notation simple, we commence by choosing the statistical 
interaction to be fermionic at first. The generalization to generic Haldane statistics 
is straightforward thereafter. First of all, we Fourier transform the TBA-equations 
in the variant (|^) 

2nlik,r)=j2^^,ik)L,ik,r). (71) 

i=i 
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For the reason that the Fourier transforms of rmjcosh^ and ei{6,r) do not exist 
separately, we introduced here the quantity ^^{0,r) = rmicosh.9 —ei{9,r). This 
should be kept in mind before transforming back to the original variables. Substi- 
tuting now the general form of ^ij{k) from (|^) into equation (^) and taking the 
inverse Fourier transformation thereafter yields 






Here we introduced the universal kernels 

h 



(72) 



^h{e) 



2 cosh f^ 



^ ^ ,^, 2/1 sin ^5 cosh §^ 

and nh,Bid) = T-d ^ 

cosh no — cos ttB 



(73) 



Neglecting the coupling constant dependent term involving VLh^sid)-, i-e. taking 
the limit B — >• ico, the identities (0) coincide precisely with the equations (7) in 
0]. This is to be expected since the latter equations describes the system which 
dynamically interacts via the minimal part of (|5(]|). Note further that the strong- 
weak duality is still preserved, i.e. Qfi,B{(^) = ^h,2-B{d)- 

We may now easily derive the generalized form of ( |72D for generic types of 
Haldane statistics. Recalling that for the derivation of (H) all the L— functions 
were parameterized by Ci{9,r) = In f 1 + e~^'(^'^M and the statistical interaction by 
9ij — ^ij " 9ijy ^^6 Fourier transformation reads 



Uk,r) = Y: [g[^q{k,r) + ^^^{k)q{k,r)/2n 



(74) 



i=i 



In the same manner as for the fermionic case we derive 



uo,r)=Y: 



h,i^,-Eih + 9ji)Q)*n, + g[^C, 



i9,r) + {Q*nh,B){e,r) . (75) 



Of course we recover (|7^) from ([7^) for g'^^ = 0. 

As a final consistency check we carry out the limit to the free theory, that is 
5^0. We have 



lim{Ci*nh,B)id,r) 



lini - 

B-.0 2 



de' 



hB cosh l^e' 



cosh /i^' - 1 + 7r25 2/2 



B 



C,{e-e',r) (76) 



/ i-'TT 

dt lim — -C: i 6 — — arsinh — , r 

B-o vr t2 + 52 ' V h 2 ' , 



C.{0,r) 



(77) 
(78) 
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In the last equality we employed the well known representation for the delta function 
tt6{x) = lini£_>oe/(e^ + x"^). This means that in the extreme weak limit the TBA- 
equations (ffSD is solved by 



ei{9, r) = rnii cosh 61 - Y^{6ij + g[-)Lj{6, r) . 



(79) 



i=i 



That (179| ) are indeed the TBA-equations may also be seen by taking this limes 
directly in (|[) and noting that \miB->o Viji.^) = ^iji^^)- In this way we have "switched 
off" the dynamical interaction and obtained a system which interacts purely by 
statistics. In particular for fermionic type of statistics {g[j = 0) or for g'^j = —6ij 
these equations describe r free bosons or r free fermions, respectively. Taking 
g'ij to be diagonal, the equations ([7^) are equivalent to a system whose S-matrix 
corresponds to the one of the Calogero- Sutherland model and whose statistical 
interaction is taken to be of fermionic type 



ijrrij 



2mi cos{7i / h) [^, from 
rrrii cosh 6' = r ^ ^ij^^j (cosh *Qh) {&■, ^) 



We finish this subsection by noting that X]j I; 
which follows 



such that the entire r-dependence and the mass spectrum in 
eliminated. We obtain 



e^{e,r) = Y. 



3=1 L 



hi^i + Y.i^ji + 9'ji)C,i) *^h- g[jCj 



1=1 



iO,r) 



^ 




(80) 


m 


and (|75|) 


may be 


{C,* 


^h,B){0/ 


r) (81) 



instead. One may solve these equations and re-introduce the r-dependence there- 
after by means of the asymptotic behaviour. This is similar to the situation in 
the next subsection, in which the equations are even further simplified to a set of 
functional equations. 



4.4 Y-systems 

An alternative analytical approach is provided by exploiting properties of what is 
often referred to as Y-systems 0. In order to derive these equations, we invoke 
once more the convolution theorem on the convolution in equation (P) and compute 
thereafter the sum of (||) at O + iir/h and 6 — iir/h minus J2j hj times (H). Using once 
more the fact that ^jIijTnj = 2mi cos{tt / h) , the terms involving rrii cancel. Em- 
ploying thereafter the identity ( ^3]) and transforming back into the original 6'— space 
we obtain, upon the introduction of the new variable Fj {9) = exp{—ei{9)), 



YA9+f]YA9-f 



l + Yii9 + f{l-B)] l + Y,l9-f{l-B 



X n 1 + Yr' (9) 
j=i V 



iji 



1 + Y, ^+f 



l + Y,[9-f)\f\[l + Yi 



-/^,\-^ 
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For fermionic statistics, i.e. g^^ = we recover the systems stated in [Q {h = 2) and 



[ TsU , when we chose the "roaming parameter" therein such that S{6) equals (pOf ). 
The strong-weak duahty is evidently still preserved. The Y-systems corresponding 
to the "minimal theories" are obtained by taking the limit B — >■ ioo, by noting 
that limg^oo Yt {0) = 0. The latter limes follows from the TBA-equations together 
with the explicit form of the TB A- kernel, such that for large 9 we have ei{9) ~ 
rTJii cosh 6'. For fermionic type of statistics the functional equations obtained in this 
limit coincide with those found in [^. As is to be expected from (^), we obtain for 
g[j = (for h = this holds in general) in the extremely weak coupling limit 

Y, {9+2m) = Yi{9), B ^ 0,2 . (82) 

To derive these periodicities directly from the functional equations is rather cum- 
bersome once the algebras are more complicated and we do not have a general case 
independent proof. This is similar to the situation in which the system involved 
is the minimal part of ( PU] ) [Q. We leave it for future investigations to exploit the 
functional equations further. 

4.5 Scaling functions 

As already stated, solving the TBA-equations and the subsequent evaluation of 
(^ leads to expressions of the scaling function c{r) with the inverse temperature r 
times a mass scale as the scaling parameter. We may carry out this task numerically. 
Alternatively, assembling our results from section 3 and the analysis of the TBA- 
kernel, we can write down a universal formula for the scaling function of the system 
in which the statistical interaction is of fermionic type and in which the dynamical 
interaction is governed by the scattering matrix of affine Toda field theories related 
to simply laced Lie algebras, up to the first leading order in the ultraviolet regime. 
From (p9D and (|6^) it follows under the assumption that the constants /3j are the 
same for all particle types 

. f n^B(2-B)(h + l)\ , , 

We will now compute step by step several scaling functions for affine Toda field 
theories in which we specify some concrete algebras. We gradually increase the 
complexity of the models and focus on different aspects. 

4.5.1 aS^^-ATFT (Sinh-Gordon) 

The Sinh-Gordon theory is one of the simplest quantum field theories in 1-1-1 di- 
mensions and therefore the ideal starting point to concretely test general statements 
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made in this context. Semi-classical investigations of a TBA-type nature concern- 



ing this model may already be found in |3l|. The model contains only one stable 
particle which does not fuse and its scattering matrix simply reads 



s'^ie) = {i}, 



tanh ^ (9 



fB 



tanh i (^ + f 5 



U) 



For a system whose dynamics is described by this S-matrix and whose statistical 
interaction is of general Haldane type, the TBA-equation (^ reads 



rm cosh 6^ + ln(l - e'^^'^'''') 



1 

2^ 



de'ip^^{e')L{e-e',r)-gL{e,r) (85) 



with 



^'^ {9) = Q,^B{e) 



4sin(7r5/2)cosh0 



6) 



cosh 29 — cos vri? 

We may convince ourselves that (^) can be obtained from (j^) and the direct com- 
putation of the kernel as well as from ( [75| ) by noting that the incidence matrix of 
course vanishes for the Ai related Dynkin diagram. 
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Figure 2: Numerical solution for the Sinh-Gordon related TBA-equations for different 
statistical interactions g and fixed values of the effective coupling B = 0.4 and r = 10~^. 

The weak coupling limit 5 — i> is read off from ([79|) 



rm cosh 9 + ln(l — e 



_ P~L(9,r)^ 



+ gL{9,r)=L{9,r) 



(87) 
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The other hines which may be computed in the standard way is the hmes r ^ 0, 
leading to the constant TBA-equation 

ln(l-e"^(°'°)) = (1-(7)L(0,0) . (88) 



This is obviously compatible with equation (^) and indicates that the extreme 
ultraviolet behaviour is independent of the dynamics and is purely governed by the 
statistics. In particular we obtain for fermionic type of statistics L(0, 0) — *■ cxo as 
a solution of (|88D , such that with (pTSD the effective central charge turns out to be 
one. In general (|88|) yields < Ceg < 1 for (7 > 1. We also observe that due to the 
fact that L(0, 0) > 0, equation ( pSf ) does not possess any physical solutions at all 
for g < 1, which means in particular that there are no solutions for the statistical 
interaction of bosonic type. 

We shall now turn to the full solution of the TBA-equation (|85|) . First of all 
we solve this equation numerically. For (7 > 1 we observe in figure 2 the typical 
behaviour, known from the analysis of the minimal part of the S-matrices 0, that 
the L-function is essentially constant between ±21n(r/2). As expected from the 
previous discussion, also the numerical analysis does not yield any solution for g < 1. 
For a more detailed discussion concerning the nature of the numerical procedure 
and in particular the question of convergence and uniqueness of the solutions we 
refer to section 5. 

We shall now compare our numerical results with our analytic approximations 
of section 3. The general formulae (|45| ) or (|48D yields for the Fourier transform of 
the TBA-kernel 

27rcosh f— - ^^^1 

$^^(A;) = ^l ' ^ - 27rg . (89) 

cosh(f) 

For the reasons mentioned in section 3.2, we will now restrict ourselves to 
fermionic type of statistics, such that 

r7(°) = and ^(2) ^ ^'^{2 - B) ^^^^ 

8 

We obtain a fairly good agreement between the numerical solution of @) for the 
function L{6,r) and the approximated analytical solution (p^), if we choose the 
constants in (|20|) to be 

/? = ln(i?(2- 5)2^+^(2-^)) and a= "" , (91) 

2(/? - ln(r/2)) 

such that 

L{e, r) = In (cos^ (aO)) - ln[5(2 - B){TTa/2f] . (92) 
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Figure 3: Numerical solution (a) versus approximated analytical solution (b) for the Sinh- 
Gordon related TBA-equation for various values of r and fixed effective coupling B = 0.4. 



Figure 3 shows that the r dependence of L{6,r) is captured very well by the 
approximated analytical solution (P3). As expected from the nature of the assump- 



tions made in the derivation, (p^) becomes relatively poor when L{6,r) is small. 
With regard to the aim of these computations, that is the evaluation of the scaling 
function c{r), precisely in this region the error is negligible as is seen from (Q). No- 
tice also that the assumptions (|16D are justified in retrospective by the numerical 
solutions. The quality of the approximation may also be seen by comparing tables 
1-3 for large values of L{9, r). 
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Figure 4: Numerical solution versus approximated analytical solution for the scaling 
function of the A\ ,A2 and Ag -affine Toda field theory related TBA-systems with 
fixed effective coupling B = 0.4. 

Finally we compute the scaling function, which acquires in this approximation 
the form 



c^G(r) 



3 7r^B{2 - B) 



8(/3 - ln(r/2))2 ' ^ ' 

Also this function is well approximated in the ultraviolet regime (r < 1) as is seen 
from figure 4 and the tables 1-3. 



r 


L(0,r)/n. 


L(0,r)/a. 


c{r)/ n. 


c(r)/a. 


10-1 


1.845 


2.135 


0.8767 


0.9281 


10-2 


2.957 


3.238 


0.9648 


0.9762 


10-3 


3.722 


3.945 


0.9845 


0.9882 


10-4 


4.286 


4.467 


0.9914 


0.9930 


10-5 


4.740 


4.880 


0.9949 


0.9954 


10-6 


5.153 


5.222 


0.9975 


0.9967 



Table 1: Numerical solution (n.) versus approximated analytical solution (a.) 
Sinh-Gordon related TBA-system with fixed effective coupling B = 0.1. 



for the 



23 



r 


L(0)/n. 


L(0)/a. 


c{r)/ n. 


c(r)/a. 


10-1 


1.478 


1.249 


0.8311 


0.8257 


10-2 


2.307 


2.220 


0.9339 


0.9340 


10-3 


2.915 


2.870 


0.9654 


0.9655 


10-4 


3.387 


3.361 


0.9789 


0.9789 


10-5 


3.771 


3.754 


0.9857 


0.9858 


io-« 


4.094 


4.082 


0.9897 


0.9897 



Table 2: Numerical solution (n.) versus approximated analytical solution (a.) for the 
Sinh-Gordon related TBA-system with fixed effective coupling B = 0.4. 



r 


L(0,r)/n. 


L(0,r)/a. 


c{t)/ n. 


c(r)/a. 


10-1 


1.335 


0.811 


0.8070 


0.7297 


10-2 


2.074 


1.782 


0.9165 


0.8977 


10-3 


2.631 


2.433 


0.9540 


0.9466 


10-4 


3.072 


2.923 


0.9709 


0.9673 


10-5 


3.435 


3.317 


0.9800 


0.9780 


10-6 


3.743 


3.646 


0.9853 


0.9841 



Table 3: Numerical solution (n.) versus approximated analytical solution (a.) 
Sinh-Gordon related TBA-system with fixed effective coupling B = 0.9. 



for the 



4.5.2 



A?- 



l(2) 



ATFT/A^^^-ATFT (Bullough-Dodd) 



In comparison to the Sinh-Gordon model the next complication arises when we 
allow the particle in the system to fuse. This case will arise when we consider the 



iW 



I (2) 



A2 and A2 -ATFT, which are known to be closely related. The former contains 
two particles which are conjugate to each other, whereas the latter contains only 
one particle which is the bound state of itself. With regard to the TBA-analysis 
these models may be treated on the same footing. The corresponding two particle 
scattering matrices read 



-,A)- 



(1) 



-,A)- 



(1) 



l(2) 



sf{d) = S2i {e) = {l}e S:i {d) = {2}, S^^^-'{d) = {l},{2},. (94) 



l(2) 



1(1) 



l(l) 



Due to the fact that S^^" {6) = S^' (0) S^i ' (9), the TBA-equation of the A^ - 
theory equals the two TBA-equations of the A2 -theory under the natural assump- 
tions that -^i(^) = L2{0) and Qij = Qji. The common TBA-kernel without the 
statistics factor is computed to 

-4v^cosh^ 



v.--(^) 



4cosh^sin((l + 5)7r/3) 



4cosh^sin(57r/3) 



2 cosh(2^) + 1 cosh(2^) - cos((l + fi)27r/3) cosh(2^) - cos(257r/3) ' 

(95) 
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We may then write the TBA-equation either in its variant (^ as 

oo 

. cosh e + ln(l - e"^(^'")) = — / d^V^^ {6') 1(6 - 9', r) - gL{e, r) (96) 



rrrii 



or in its universal form ([75|) as 

aO, r) = [(e - (1 + 9')C) * fig] (^, r) + g'L{B, r) + (C, * fi,,^) (e, r) . (97) 

We used the abbreviations g = gn+gu and g' = g'li+g'u- The Fourier transformed 
TBA-kernel is computed from (^Sj) or (^Sf) 



_.(2) -4(1) ~/l(l) -4(1) -4(1) 

47r fsinh f^) + sinh f^)) fcosh (^ - ^) - i) 

= ^ ^^ ^ ' /,\ ^ ^^ ^ - 27r(7 . (99) 

sinh(7rA;) 

We restrict now to fermionic type of statistics and chose the constants in our analytic 



solution (pOD in the same way as for the Sinh-Gordon model, i.e. 



/3 = /3^ = /52 = ln(5(2- 5)2^+^(2-1?)) ^^^ a = a^ 



^2 



2(/?-ln(r/2)) ' 

(100) 
such that the approximated function for L{6,r) equals (|92[) . The scaling function 
up to the first leading order (p9D becomes now 



Solving (p6| ) numerically we obtain qualitatively the same kind of agreement with 
the approximated analytical solution as for the Sinh-Gordon model, see figure 4. 

4.5.3 A^^^-ATFT 

The simplest model which involves two inequivalent TBA-equations coupled to each 
other is the A^ -afiine Toda field theory. In our conventions the two-particle scat- 
tering matrices for this theory read 

Sf {0) = Sf {e) = {l},, Sf (0) = {3}, , (102) 

Sf (0) = sf (0) = {2}, , Sti' (0) = {1}, {3}, . (103) 



Particle 1 is the anti-particle of 3 and particle 2 is self-conjugate. The masses are 
mi = m^ = m/\/2 and m2 = m. In the direct channel the fusings 11 ^ 2 and 
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33 — > 2 are possible. Under the assumption that ei{9,r) = 63(6, r), the TBA- 
equations in the variant (§) now read 



^(2) aW 



ei{9,r) = rm/ V^ cosh.9 — {lp22 * ^1 + fii * ^2 

-ig[,+g[,)C,ie,r)-g[,C2{e,r), (104) 

/ ^(2) ^(2) \ 

62(9, r) = rm coshO — {ip22 * C2 + 2ipi.f * Ci\ {9,r) 

-W12 + 9'2^)'^i{0, r) - g'22C2{e, r) , (105) 

where the kernels are 

Af ^ MB7r/4:) cosjBn/A) 2_ 

^^^ ^ ' cosh^^-cos2(57r/4) cosh^^ + sin2(S7r/4) cosh^ ^ ^ 

4'Vm o ^^^ 2sin((i?+l)7r/4) ^2_\ 

^- (') = '^"^^H cosh(2^)+sin(i?^/2) "^^i^j • ^ °'^ 

.(2) .(2) -(2) .(2) .(2) 

Note that ipfi {9) = (^23 (^) and (^22 (^) = ^n {0) + (^fl (^)- Ahernatively 
we may write down the universal form of the TBA-equations (0), which becomes 
particularly simple for the fermionic statistics 

UO,r) = {^2-L2)*n,{9,r) + {L,*n,,B){9,r) (108) 

^2{9,r) = 2{i,-L^)*n^{9,r) + {L2*n^,B){9,r) . (109) 



We compute from (^51) or (^Sj) the Fourier transformed TBA-kernels 



~Af\ . _ (1 + 2cosh(7rfc/4))(cosh((E - l)7rA:/4) - cosh7rfc/4) 
^22 ^^^ " ^"^ cosh(A;7r/4) + cosh(3A:7r/4) ^^^^^ 

42) ^ ^^cosh(../4)cosh((E-l)../4)-l _ 

cosh(/c7r/2) 

From (|60|) or from the explicit expression for /i^!^^ corresponding to each block in 
( |102| ) and ( |103D , we obtain for fermionic statistics 




„<»' = and ,g' = "'''''i' "" I 2 4 2 I . (112) 



Therefore we have 



(2) (2) 37r^i?(2 - B) (2) Ti^B{2 - B) 
m =m = ^ and r]\ = • (113) 
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/r)\ /rf\ fr\\ 

Computing the sum ^1+^2+^3 confirms our general formula (|64D. Finally we 
obtain for the approximated scaling function 

Once more we compare this analytical expression with the numerical solution of 
the two coupled TBA-equations (|104|) and ( |105|) . Choosing the constants in the 
same way as for the Sinh-Gordon model, figure 4 demonstrates the same kind of 
qualitative agreement as we observed for the previous models. For all models we 
observe that at about r = 0.1 the analytical expressions for c{r) start to decrease 
more rapidly than the numerical solution. 

4.5.4 (G^^^ ^ Df^)-ATFT 

This theory is the simplest example for a model in which the masses as a function 
of the coupling constant flow from the classical mass spectrum of one Lagrangian to 
its dual in the Lie algebraic sense. The theory contains two self-conjugate particles 
whose masses were conjectured ||S^, ^ (supported by numerical investigations P^ ) 



to rrii = msm{7T/H) and m2 = msm{27r/H). The "floating Coxeter number" is 
taken to be iJ = 6 + 3B with < i? < 2. The fusing processes 11 — > 1 + 2, 
22 —>■ 1 + 2 + 2 and 12 -^ 1 + 2 are possible. In our conventions the related 



scattering matrices |32, 26 1 read 



^i7^ (9) = [1, 0]e[H - 1, 0]e[H/2, 1/2],, Sf/^ (9) = [H/3, l]e[2H/3, l]e, (115) 
SS"" (9) = [H/3 - 1, l]e[H/3 + 1, l]e[2H/3 - 1, l]e[2H/3 + 1, 1], . (116) 

Besides the loss of duality these scattering matrices exhibit a further difference in 
comparison with the simply laced case, which has a bearing on the TBA-analysis. 
In the standard prescription the symmetry of the Bethe wave function is derived by 
exploiting the behaviour of the scattering matrices at ^ = 0. As already mentioned 
for the simply laced case we always have Sij{0) = (— 1)^*^, whereas now we observe 

'S'22 (0) = "'S'j^i (0) = 1- Assuming that the particles described are bosons, 
we should choose, according to the arguments of 0, (711 = 1 and (722 = 0. In 
this case, however, the TBA-analysis does not produce any physical solution. The 
TBA-equations in the variant (H) become 

6i{9,r) = rm sm{TT / H) cosh9 — ((f II * CA {9,r) — ((fi2 * ^2) i9,r) 

-g[,C,i9,r)-g[,C2i9,r), (117) 

£2(9, r) = rm sm{27c/H) cosh 9 - (^^^ * Ci) {9, r) - (^^^ * £2) {9, r) 

-g'^,C,{9,r)-g',,C2i9,r) . (118) 
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The kernels and their Fourier transformed versions read 



and 



^Tio) = 


UJi{e) + UJh-i{6) + UJh/2,1/2{0) 


(119) 


^T{0) = 


V'2l(^) = ^H/3,l(^) + C^2H/3,l(^) 


(120) 


^T{0) = 


^H/3-l,l(^) + ^H/3+l,l(^) + ^2H/3-l,l(^) + ^2H/?,+l,l{d) 


(121) 


^Tik) 


cosh((l + B)'Kk/H) - cosh((l - 2B)7Tk/H) 


(122) 


1-2 cosh(7rA;/3) 


^Tik) 


^12^ (k) + 47r cosh((l - B)nk/H) 
2cos\i{Txk/H) 


(123) 


^T{k) 


= 2 cosh{n k/H)^f^^{k) + 2n . 


(124) 



From this or our general formulae of section 4.2. we obtain 

(2) 7r'^B{2 - B) (2) (2) (2) „ (2) (2) ^ (2) „ (2) 

^11 = jp ' ^12 = -^11 ' ^22 = -2^11 , ^1=0, r]\' = -?>ri\{ . 

(125) 
The facts that rj\ vanishes and rj2 becomes negative make it impossible to 
use our approximated analytical solution (pOf). From a numerical point of view not 
much has changed in comparison with the previous models and we can still solve 
the related TBA-equations in the standard way. The related L-functions look qual- 
itatively the same as in the previous cases and the scaling functions for some values 
of the effective coupling are depicted in figure 5 together with the same function for 
the A2 -affine Toda field theory. We observe that for fixed value of r, the scaling 
function for the (G2 > -D4 )-affine Toda field theory is a monotonically increasing 
function of B. Under the same circumstances the scaling functions related to theo- 
ries in which the underlying algebras are simply laced are monotonically decreasing 
at first up to the self-dual point B = 1 and monotonically increase thereafter. 
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Figure 5: Scaling functions for the A>2' and (G2 ,1)4 '')-affine Toda field theory related 
TBA-systems for various values of the effective coupling. 



5 Existence and Uniqueness Properties 

In this section we are going to investigate the existence and uniqueness properties 
of the solutions of the TBA equations. Our main physical motivation for this con- 
siderations is to clarify whether it is possible to obtain different effective central 
charges for a fixed dynamical and statistical interaction due to the existence of sev- 
eral different solutions of the TBA equation. As a side product we obtain useful 
estimates on the error and the rate of convergence of the applied numerical proce- 
dure. Precise estimates of this kind were not obtained previously in this context 
and convergence is simply presumed. The method we employ is the contraction 
principle (or Banach fixed point theorem), see e.g. [^. For the Yang- Yang equa- 
tion in the case of the non-relativistic one-dimensional Bose gas for fermionic type 
of statistics the uniqueness question was already addressed by Yang and Yang |l[, 
albeit with a different method. 

In order to keep the notation simple we commence our discussion for a system 
with one particle only and a statistics of fermionic type. Thereafter we discuss the 
straightforward generalization. The standard way to solve integral equations of the 
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type (ID consists in discretising the original function and a subsequent iteration, 
usually numerically. Normalizing the mass to one, this means we consider (^ as 

en+i{e,r) := r coshe - cp *\n{l + e-'"^'^''^^) (126) 

and perform the iteration starting with 6Q{6,r) = rcosh^. The exact solution is 
then thought to be the limes lim„^oo£^n- However, a priori it is not clear whether 
this limes exists at all and how it depends on the initial value Eq. In particular, 
different initial values might lead to different solutions. 

The natural mathematical setup for this type of problem is to re-write the TBA- 
equation as 

iAOie,r) :=(^*ln(l + e«(^'^)-'-™^^^) =e(^,r), (127) 

and treat it as a fixed point problem for the operator A. 

In order to give meaning to the limes lim„^oo ^n we have to specify a norm. Of 
course it is natural to assume that C, as function of 6 is measurable, continuous and 
essentially bounded on the whole real line. The latter assumption is supported by 



all known numerical results. Furthermore it follows from ( |127D , together with the 



explicit form of (p, that possible solutions ^ vanish at infinity. This means possible 



solutions of (|127D constitute a Banach space with respect to the norm 



, = esssup 1/(^)1, (128) 

i.e. Loo(lR). In principle we are now in a position to apply the Banach fixed point 
theoremlH], which states the following: 

Let D C Loo be a nonempty set in a Banach space and let A be an operator 
which maps D g-contractively into itself, i.e. for all f,g & D and some fixed q, 
< g< 1 

\\Aif)-Aig)\\^<q\\f-g\\^. (129) 

Then the following statements holds: 

i) There exists a unique fixed point ^ in D, i.e. equation ( |127| ) has exactly one 
solution. 

ii) The sequence constructed in ( [126| ) by iteration converges to the solution of 



iii) The error of the iterative procedure may be estimated by 

||^-^„||oo<^ -ll^i-^olloo and ^-^„+i <^ ^n+l-^r. 

L — y oo j^ — Q 



" One may of course apply different types of fixed point theorems exploiting different properties 
of the operator A. For instance if A is shown to be compact one can employ the Leray-Schauder 
fixed point theorem. In the second refence of B it is claimed that the problem at hand was treated 
in this manner, albeit a proof was not provided. 
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iv) The rate of convergence is determined by 

^-^n+l\\^<QU-U\oo- 

In order to be able to apply the theorem we first have to choose a suitable 
set in the Banach space. We choose some q G [0, 1) such that e~^' < q and take 
D to be the convex[2] set Dqr := {/ : ll/lloo < In yz~ + r\. We may now apply 
the following estimate for the convolution operator (f* (which is a special case of 
Young's inequality) 

||<^*/||oo<||<^||l||/||oo. (130) 



For the concrete one particle models at hand, the Sinh-Gordon- and the Bullough- 
Dodd-model, we have ||v9||i := J ^\f{0)\ = 1. For the operator L we have the 
estimate 



L{f) < In [1 + exp 



- r)l < In < In — ^ + 

'^ - 1-a - 1-0 



r. 



The last inequality follows from our special choice of q. Thus, A maps Dg,. into 
itself. 

In the final step we show that the contraction property ( |129| ) is fulfilled on Dgj. . 
It suffices to prove this for the map L, because of ( |130| ) and the fact that \\(p\\i = 1. 
We have 



ll^(/)-^(^)l|c 



}dt^L{g + t{f~g)) 
at 



Jdt 




if -9) 






1 + exp(-5f -t{f - g)+r cc 
1 


)sh9) 

\f- 


oo 
^loo 

\f-9 


0<t<l 


1 + exp{-g - t{f -g) + r) 
1 


0<t<l 


l + expi-\\g-t{f-9)\^ 


+ r) 



< q\\f-9\L ■ 

In the last inequality we used the fact that Dgr is a convex set. 

We may now safely apply the fixed point theorem. First of all we conclude from 
i) and ii) that a solution of ( |127| ) not only exists, but it is also unique. In addition 
we can use iii)and iv) as a criterium for error estimates. From our special choice 
of the closed set Dg^r one sees that the rate of convergence depends crucially on 
the parameter r, the smaller r the greater q is, whence the sequence (^„) converges 
slower. 



-For f,gG Dg^r also t/ + (1 - t)g G Dg^r, with < t < 1. 
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One could be inathematically more pedantic at this point and think about dif- 
ferent requirements on the function ^. For instance one might allow functions which 
are not bounded (we do not know of any example except when r —>■ 0) and then 
pursue similar arguments as before on Lp rather than Loo- 

The generalization of the presented arguments to a situation involving I dif- 
ferent types of particles and general Haldane statistics may be carried out in a 
straightforward manner. 

6 Conclusions 

Clearly an unsatisfactory feature of our approximated analytical expression for the 
scaling function (|83| ) is that, at the moment, we do not have any additional con- 
straint at hand which allows to determine the constant f3. Admittedly the choice 
(pip enters our analysis in a rather ad hoc way. At the moment the rational behind 
our choice is that we may compare it with the semi-classical results in the spirit of 



T4 and thereafter restore the strong-weak duality. In addition it is supported by 
our numerical results. Surely this is by no means unique and it is highly desirable 
to eliminate this ambiguity. 

Nonetheless, different choices of the constant will not change the overall be- 
haviour and our expression is clearly in conflict with the results of Cassi and Destri 

T6[| , who found /(r) = —3r'^ml/{7r^[i) for the function in equation (p. The con- 



vention therein are that 1 labels the particle type with smallest mass and ipn results 
from the power series ^PuiO) = — Z^fcli V^ii exp(— fc \6\). On the side of the TBA- 
analysis the origin of this discrepancy may be tracked down easily. A behaviour of 



the type quoted in [|T6| was derived before in [|T^, ^ for afiine Toda field theories 
related to the minimal part of the scattering matrix (pOl). This result was then 
extrapolated by Cassi and Destri in the way that ipii was taken to be the coeffi- 
cient of the power series related to the full coupling constant dependent scattering 
matrix. However, the derivations in |1T2| , ^ rely heavily on the assumption that the 



function L{9) is constant in the region — ln(2/r) ^ ^ ^ ln(2/r). It is essentially 
this property which is lost for the full theory (as our numerical results demonstrate), 
such that the arguments of Cassi and Destri become faulty. It would be interesting 
to settle the question also on the perturbative side and bring our results into agree- 
ment with perturbation theory. Our findings will surely have consequences for the 
subtraction scheme used in such considerations. 

There exist of course other methods which allow to extract the ultraviolet central 
charge from a massive integrable model. The c-theorem |]3B[ has turned out to be 
very efficient in this context and a direct comparison is very suggestive. As the 
main input, the c-theorem requires the n-particle form factors (the 2-particle form 
factor is usually sufficient) related to the trace of the energy-momentum tensor. 
Unlike in the situation of conformal invariance, in which the trace vanishes, this 
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tensor is not unique and acquires some scaling behaviour for massive theories. As 



observed in [^, the consequence of this fact is that the c-theorem produces a whole 
ray of central charges greater than the rank of the underlying Lie algebra g. The 
natural question which arises is to identify the origin of this ambiguity inside the 
TBA-approach. One might assume that it results from several possible solutions 
of the TBA-equation. However, this possibility is definitely ruled out as follows 
from our investigations in section 5. Since the statistical interaction g enters the 
analysis as a further parameter, this could provide a further possible mechanism 
which produces the observed values. As our investigations demonstrate, however, 
also this possibility can be excluded for certain, since different choices for g only 
produce central charges smaller than the rank of g. In the light of this results we 
conjecture that the responsible mechanism is to include a non-vanishing chemical 
potential in the way as was already indicated by Yang and Yang [|^ . To settle this 
question precisely requires more detailed investigations [^ . 



There are several further questions which should be addressed in order to com- 
plete the picture. It will certainly be interesting to obtain also the higher terms 
beyond the first leading order in (|l]) and to exploit further the Y-systems of section 
4.4. as an alternative analytical approach. There exist interesting links between 
these systems and spectral functions |3^, such that one can expect more exact and 



universal results to follow. 
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